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We present the first thermoelectric measurements on pentacene field effect transistors. We report 
high values of the Seebeck coefficient at room temperature between 240 and 500/iV/K depending 
on the dielectric surface treatment. These values are independent of the thickness of the channel 
and of the applied gate voltage. Two contributions to the Seebeck coefficient are clearly identified: 
the expected contribution that is dependent on the position of the transport level and reflects 
the activated character of carrier generation, and an unexpected intrinsic contribution of 265 ± 
40 A*V/K that is independent of the temperature and the treatment of the oxide surface. This value 
corresponds to an unusually large lattice vibrational entropy of 3/cb per carrier. We demonstrate 
that this intrinsic vibrational entropy arises from lattice hardening induced by the presence of the 
charge-carrier. Our results provide direct evidence of the importance of electronic polarization effects 
on charge transport in organic molecular semiconductors. 

PACS numbers: 72.80 Le - 72.20 Pa - 33.15 Kr 



I. INTRODUCTION 

Due to the narrow-band nature of their electronic 
states, charge carriers in organic molecular semicon- 
ductors are prone to perturbations by electronic and 
lattice excitations^ Despite decades of research devoted 
to understanding the implications of these interactions 
to charge transport, a clear consensus is yet to emerge. 
The current wide interest in these materials for potential 
optoelectronic applications (see the review articles in 
Ref. 1 exacerbates the need for gaining insight into 
this problem. Recent papers have emphasized the 
role of transfer-integral fluctuations^ or polarization 
fluctuations^ in the transport process, particularly in 
the context of experimental observations of the decline 
of the carrier mobility with increasing gate-dielectric 
permittivity seen in both single-crystal and thin-film 
organic transistor devices*^ 

Thermoelectric measurements, because they provide 
direct access to the heat transported by the carrier, can 
unravel the origin of charge transport. Surprisingly, 
to our knowledge, Seebeck coefficients have never been 
measured in the crystalline organic semiconductors such 
as pentacene, rubrene and related compounds of interest 
here. 

In this article, we present what, to our knowledge, are 
the first results of thermopower measurements on pen- 
tacene thin transistors. We find that a temperature- 
independent, intrinsic contribution to the Seebeck coeffi- 
cient of 265 ± 40 /uV/K is present in all devices regardless 
of the gate-insulator used. This value corresponds to a 
very large lattice vibrational entropy of about 3fcs per 
carrier. We shall show that this result arises from the 
creation of an intermolecular phonon cloud associated 
with local changes in the Van-der-Waals interactions be- 
tween molecules that accompanies the electronic polar- 



ization induced by the carrier in the channel of organic 
field-effect transistors (OTFTs). 



II. EXPERIMENTAL 

A. Fabrication procedure 

For this work, top-contact thin (5 nm, 10 nm and 
100 nm thick active layer) pentacene transistors were 
fabricated by standard high-vacuum deposition tech- 
niques as we reported previously!^ Both thermal SiC>2 
(200 nm thick) on heavily-boron-doped silicon wafer and 
single-crystal AI2O3 (sapphire) were used as the device 
gate-dielectric. 

We also explored possible intcrfacial effects associated 
with the surface treatment of the oxide layer under 
various plasma exposure conditions and by vacuum 
deposition of a self-assembled monolayer (SAM)^ prior 
to pentacene growth. Plasma treatment with argon was 
done at a radio frequency of 10 W for 5 minutes at a 
constant gas pressure of 0.1 mbar. The oxygen plasma 
treatment was done under the same conditions except 
the exposure time was set to 0.25, 2, 3 and 10 minutes. 

To form large pentacene grains within the conducting 
channels of the transistors, a substrate temperature of 
328 K and a deposition rate of 1.2 nm/min were main- 
tained during growth^ The top source and drain gold 
contacts were evaporated through a shadow mask at the 
temperature of liquid nitrogen. The channels were fab- 
ricated with a width of 6 mm and a length of cither 
50 or 100 /Urn. The four-probe contact geometry con- 
sisted of four 6 mm- wide gold strips spaced 600 [ira apart. 
The morphology of the evaporated films was studied by 
atomic force microscopy (NT-MDT Solver Pro in tapping 
mode). 
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FIG. 1: Schematic diagram of thermoelectric measurement 
apparatus. The sample was fixed by a Apezion® high vacuum 
grease to the copper heating elements to ensure the thermal 
contact. The thermocouples were cemented with Stycast® on 
the outer gold pads. The electromotive force AV was mea- 
sured on the inner gold pads. The gold wires were glued with 
standard colloidal silver on the gold contacts. Thermocou- 
ples Ti, T% and T3 were used to measure the temperature of 
the heating element and of the copper mass (cryostat). The 
cryostat was electrically insulated and normally grounded en- 
suring Vg = V. 



the activation energy of the thin-film conductivity $ a . 
The I — V characteristics of our pentacene thin-film 
transistors yielded field-effect mobilities in the range of 
0.1 to 0.4cm 2 /Vs as shown in Table [I] These mobility 
are deduced from the linear region of the I — V charac- 
teristics. The activation energies are calculated from the 
slope of the associated Arrhenius plot. A typical result 
is given and analyzed in Figure [5] 
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1 


2 180 sec 


10 


0.22 


104 


95 
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Ar 300 sec 


10 


0.32 


80 


32 
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2 300 sec 


5 


0.33 




58 


4 


O2 300 sec 


100 


0.2 


84 




9 


O2 120 sec 
COOH-3acene 


10 


0.11 


161 


107 



TABLE I: The field-effect mobility \i and the activation ener- 
gies of the carrier conductivity S a and the mobility of our 
pentacene thin-film transistors. 



B. Characterization 

Dark current versus voltage (I — V) characteris- 
tics were obtained under ambient conditions using a 
Kcithlcy 4200 semiconductor characterization system. 
The thermoelectric, four-probe and activated mobility 
measurements were carried-out under high vacuum of 
1 x 10 -6 mbar with a home made vacuum chamber. 

Thermoelectric power was measured by fixing the sam- 
ple to two copper heating blocks (see Figure [T|) . The 
temperature was measured with a differential NiCr/CuNi 
thermocouple adhered to the gold contacts with a ther- 
mally stable resin. A temperature gradient of 5 K along 
the sample was imposed for the duration of 1 minute. 
The Secbcck coefficient was derived from the electromo- 
tive force AV divided by the applied temperature gra- 
dient AT. To ensure that there were no artifacts due 
to the experimental setup, the direction of the temper- 
ature gradient was reversed and the measurement was 
repeated. 



III. RESULTS 

A. Field-effect mobility and activation energy 

The electrical characteristics of an OFET provide 
insight into the transport mechanisms via the field-effect 
mobility fi, the activation energy of the mobility S^, and 
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FIG. 2: Transfer characteristics of transistor 9 generated from 
the linear regime of the I—V data in the inset. The field-effect 
mobility is /i — 0.11 cm 2 /V.s. 

The thin-film conductivities a were measured on the 
same substrate (from the same batch) by means of 
temperature-dependent four-probe technique. The acti- 
vation energy is then derived from the slop of an Ar- 
rhenius plot as shown in Figure [3l Both, the field- 
effect mobility and the channel conductivity exhibit typ- 
ical temperature-activated behavior. This suggest that 
charge transport in organic thin-film transistors occurs 
via hopping between localized states. The significance 
of our channel conductivity results is discussed in more 
detail in Section IIII CI 
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FIG. 3: On transistor 9, the field effect mobility is activated 
with the energy = 107 meV. The inset presents the four 
probe measurement of the temperature dependence of the 
channel conductivity. It is activated with S a = 161 meV . 
The differences between the Fermi level Sf and the transport 
level Sv is then £f — £v = 56meV. 



B. The thermoelectric power 

For the thermoelectric experiment, the quantity of 
interest is the Seebeck coefficient (a) which measures 
the entropy transferred to the carrier in the channel per 
unit charge. It is defined as (AV/AT), where AV^ is 
the voltage induced by an applied temperature gradient 
AT across the sample. We can neglect the contribution 
of the contacts to the thermoelectric power because the 
measured values are two orders of magnitude higher than 
that of gold. Since we measured the electromotive force 
AV in open circuit, no electrical current was able to 
pass through the thin-film. Therefore the measurement 
is independent of the geometrical shape of the contacts 
and of the contact resistance and grain boundaries do 
not have a significant influence on the results. 

The Seebeck coefficients that we obtained from our 
devices are presented in Figure [H We find surprisingly 
high values of a ranging from 240 to 500 /LtV/K at room 
temperature that depend on the nature of the device 
substrate. 

Although the results in Figure |4] show a strong depen- 
dence of the thermoelectric power on the treatment of the 
oxide surface, we do not observe a dependence on the con- 
tact resistance or the morphology of the pentacene layer. 
Table [TT] gives the measured Seebeck coefficients of pen- 
tacene on oxygen plasma activated Si/SiC>2 substrates at 
room temperature for pentacene layer thicknesses of 5, 
10, and 100 nm. It has been shown that the contact re- 
sistance of top-contact pentacene field-effect transistors 
increases with increasing active layer thickness^ Since 
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FIG. 4: Seebeck coefficients measured on pentacene films de- 
posited on the different substrates described in Table IIVI 



we observe little or no dependence of the thermoelectric 
power on the device thickness, we conclude that our mea- 
surements are independent of the contact resistance. 



Sample 


3 


1 


4 


Thickness (nm) 
Seebeck coefficient (/iV/K) 


5 
245 


10 

240 


100 
249 



TABLE II: Thermopower measurements at room temperature 
for different sample thicknesses. 

The thermoelectric power is also independent of the 
morphology of the pentacene surface prior to gold con- 
tact deposition. Figure O shows AFM images of the pen- 
tacene surface in the device channel for various film thick- 
nesses. The pentacene films shown in panels a, b, and 
c were evaporated on top of Si/SiC>2 substrates treated 
with oxygen plasma and are 5, 10, and 100 nm thick, 
respectively. Panel d shows a 10 nm-thick pentacene 
layer on top of a Si/Si02/Al203 substrate. The sub- 
strate was treated with oxygen plasma and then modi- 
fied by a COOH- Anthracene SAM. As the film thickness 
increases, growth is less homogeneous and the morphol- 
ogy of the pentacene surface is more three-dimensional. 
However, given the lack of thickness dependence we ob- 
serve in the Seebeck coefficient (Table HI)) , these changes 
in morphology do not affect the thermoelectric power. 
This conclusion and the lack of observed contact resis- 
tance are consistent with previous results that show that 
charge transport occurs in the first few molecular mono- 
layers adjacent to the gate dielectric^ 

It is also important to point out that the thermoelec- 
tric power does not depend on the gate bias strength. We 
measured the thermoelectric power of an argon plasma 
treated (transistor 2) and an oxygen plasma treated de- 
vice (transistor 3) while applying a gate potential. The 
obtained Seebeck coefficients reported in Table IIIII and 
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FIG. 5: The topography of the pentacene films with various 
film thicknesses. Pentacene films on oxygen plasma activated 
Si/Si02 substrates are shown with thicknesses of: a) 5 nm, 
b) 10 nm, and c) 100 nm. In panel d, the gate dielectric was 
modified by a SAM (COOH- Anthracene) prior to pentacene 
evaporation. The pentacene layer thickness in panel d is 10 
nm. All scans have dimensions of 10 x 10 /im 2 . 



Figure [6] do not exhibit a significant dependence on the 
applied gate potential. This confirms that the Seebeck 
coefficient in our devices is essentially a single carrier 
property.— This is true in general except in the case of 
highly correlated electron systems, where it may depend 
slightly on the carrier densityJ^ 





a (/xV/K) at room temperature 


Sample 


V G = -40 V 


Vg = V 


V G = +30 V 


2 




419 


406 


3 


221 


245 


237 



TABLE III: Dependence of the Seebeck coefficient a(/iV/K) 
on the the applied gate bias Vg- 



C. Discussion 

We shall now examine the factors that contribute to 
such high values of a. 

Following David Emm,— we recall that the main con- 
tributions to the Seebeck coefficient may be written as 
the sum of a formation contribution a f and a migration 
contribution a rn . 
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FIG. 6: The Seebeck coefficient as a function of tempera- 
ture T and applied gate voltage Vg- Transistor 2 was treated 
with argon plasma and transistor 3 was treated with oxygen 
plasma. 



Here, qaf is the change in the system's entropy in- 
duced by adding a charge carrier and is independent of 
the mechanism by which a carrier is transported through 
a material. On the other hand, a m comes from the en- 
tropy change associated with the net energy transferred 
in moving a carrier. This contribution, which depends 
on the charge transport mechanism, was calculated by 
Eminii for the case of hopping transport in disordered 
systems. In our case, a m is small when compared to a/. 
This is because of the nearly equal probability for up- 
wards hops (involving energy gain) and downwards hops 
(involving energy loss) along the carrier hopping trajec- 
tory, leading to very small values of a m of the order of a 
few ten of /xV/K. 14 i 15 If we assume further that spin con- 
tributions can also be neglected because carrier formation 
in organic semiconductors does not imply spin reorgani- 
zations, we arc finally left with two contributions that we 
can write as, 



a f 



( mixing 

q 



AS vi 



vibration ) 



(2) 



A5„ 



reflects the 



The configuration entropy, "^mixing 
changes in the distribution of the carriers between the 
electronic states. The vibrational entropy, AS v ibration, 
emanates from the modification of the vibrational modes 
due to lattice softening or hardening induced by the 
presence of a charge carrier. 

For non-degenerate hole systems, AS^ixing is well- 
know n 11 ! 16 and can be approximated by, 



a f 



(1) 



AiSVnixirier — 



T 



Sp — Sy 

T 



(3) 



5 



where Sp represents the position of the Fermi level, S 
is the hole energy and Sy is the transport energy. The 
bracket corresponds to thermal average over all the elec- 
tronic states which contribute to conduction. We can es- 
timate the value of Admixing by extracting AS = Sp—Sy 
from four-probe measurement of the temperature depen- 
dence of the channel conductivity. A typical result is 
shown in Figure [3] The quasi-two dimensional chan- 
nel conductivity portrays an activated behavior arising 
from the thermal excitation of the charge carrier from 
the transport energy to the Fermi level and from car- 
rier hopping between localized states. We modelled this 
previously^ as, 



-AS\ 

a = |e|po cxp \ -j—^- |pi exp 



|e|poMo exp 



(4) 



where S a = AS + S^ and S^ is the activation energy as- 
sociated with the field-dependent mobility. 
From this analysis, we determine the vibrational con- 
tribution to the measured Seebeck coefficients shown in 
Figure[7]as a vib rationai = Aa + a m i xin g. Table [TV] summa- 
rizes the results of our analysis. As displayed in Figure 
we find an unusually large vibrational contribution to a 
of about 265 ± 40 fiV/K that is independent of the type 
of gate-dielectric and the different surface treatment per- 
formed on the substrates. 
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FIG. 7: The temperature dependence of the mixing and vi- 
brational contributions to the Seebeck coefficient for sample 
9. 
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FIG. 8: The contribution to the Seebeck coefficient due to 
interaction of the charge carrier with its surrounding phonons 
is about 265 ± 40/iV/K and independent of the nature and 
the kind of surface treatment. 





Substrate 


Surface 


5acene 


Q 






treatment 


H 


(AiV/K) 


1 


Si/Si0 2 


2 180 sec 


10 


240 


2 


Si/Si0 2 


Ar 300 sec 


10 


419 


3 


Si/Si0 2 


2 300 sec 


5 


245 


4 


Si/Si0 2 


2 300 sec 


100 


249 


5 


A1 2 3 


2 15 sec 


10 


284 


6 


A1 2 3 


2 180 sec 


10 


283 


7 


A1 2 3 


2 600 sec 


10 


306 


8 


A1 2 3 


Ar 300 sec 


10 


429 


9 


Si/SiOa/AlaOa 


2 120 sec 
COOH-3acene 


10 


460 








Admixing 


Difference 




(meV) 


(meV) 


(MV/K) 


(AiV/K) 


1 


95 


9 


32 


251 


2 


32 


38 


172 


274 


9 


107 


56 


191 


270 



TABLE IV: Summary of the achieved thermoelectric power 
upon different substrate and surface treatments. 

"at 300 K 



When phonons with energies lower than fee are involved, 
the vibrational entropy (which in fact corresponds to the 
specific heat) follows the law of Dulong and Petit, ft 
leads to a temperature independent contribution^ given 



by, 



IV. MODEL 

We now present a model calculation that accounts for 
this very large intrinsic value of avibrationai in OTFTs. 



AS, 



vibration 



Av 



(5) 



where Nj represents the number of softened (or hard- 
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ened) modes, and Vj their frequency. 

Equation (J5J) relates the thermopower vibrational con- 
tribution to the local changes in the phonon frequencies 
induced by a charge. This contribution is exception- 
ally large for the highly polarizable organic molecular 
semiconductors^ as the nature of the chemical bonding 
is changed by the presence of a charge carrier from pure 
Van-der-Waals to charge-dipole interaction. More pre- 
cisely, if the carrier is localized on a molecule, the poten- 
tial energy in which the adjacent molecules move, is the 
sum of the potential energy arising from direct interac- 
tions between neutral molecules and the carrier electronic 
energy, itself a function of the molecule positions fl . It 
is given by, 

i 1 

where di is the dipole induced on the i th molecule by 
the carrier at the origin^ Their contribution to the po- 
larization energy § v (fi) decreases as the fourth power of 
their distance from the carrier. Then, we evaluate the 
local change in the phonon frequencies in the spirit of 
a pair-potential model. We start with a Lennard-Jones 
potential and a Van-der-Waals attractive potential for a 
pair of neutral molecules and we add the contribution of 
the carrier electronic energy^ 



For a pair of pentacene molecules, these parameters 
are estimated as follows: 7/47T£o = 40 A , ro5 A, M = 
23 x 1CT 26 kg and lu = 2itv ~ 10 13 s _1 , so that 8vjv = 
(l/2)(5k/k) = 0.32. Since the frequency hardening is 
proportional to r^~ 6 , the strongest contribution to the 
thermopower comes from the four next-nearest neighbors 
of the ionized molecule in the pentacene crystal. In our 
case, by considering two optical modes in the plane, the 
estimated thermopower is 



^vibration ~ 8 — ~ 220 /iV/K (8) 

q v 

This calculated value agrees remarkably well with the 
measured values reported in Figure [5J 



V. CONCLUSION 

The above results strongly support the conclusion 
that an electronic polaron is formed in the quasi-two 
dimensional channel of organic field-effect transistors. 
In the pentacene thin-films presented here, the charge 
carriers are localized due to the presence of intcrfa- 
cial disorder. Thus, the electronic polaron is quasi-static. 



fir) 



-Muj 2 rl 
~72 



7 



2 (47T£ ) 2 r4 V r 



(6) 



where M is the reduced mass of the molecule, tq the 
equilibrium distance between neutral molecules, 7 the 
electronic polarizability of the molecule, q is the charge of 
the carrier. The shift in the carrier equilibrium position 
from r — ro to r = Tq is obtained from the solution of 
(dy/dr) r = r * = 0. The result is, 



In high mobility transistors based on pentacene sin- 
gle crystal s 6 ' 18 the electronic polaron is delocalized as 
Bloch waves.— As further corroboration of our results, we 
predict the vibrational entropy to be weaker in single- 
crystals than in thin-films. Understanding the dynamics 
of the electronic polaron motion can lead to quantita- 
tive assessment of the observed charge transport in or- 
ganic molecular semiconductors and the optimization of 
the performance of optoelectronic devices based on them. 
We have recently proposed a new theory to predict the 
mobility of this electronic polaron 



= 1 



12 7 g 2 



(47re ) 2 ^ Mw 2 r 2 



1/6 



Mui 2 



32 79 2 



(47T£o) 2 ^ Mw 2 r 2 



(7) 



The stiffness parameter is modified from k = Mlo 2 to 
k* = k + Sk= {d 2 r/dr 2 ) r=r » and we find 
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